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Abstract 



The charged K^^ decay, tt^tt^ I'^v^ is studied in the frame- 

work of chiral perturbation theory based on the effective Lagrangian 
including mesons, photons, and leptons. We give analytic expressions 
for the two vectorial form factors, / and g, calculated at one-loop level 
in the presence of Isospin breaking effects. These expressions may then 
be used to disentangle the Isospin breaking part from the measured 
form factors and hence improve the accuracy in the determination of 
TTTT scattering parameters from Kf^^ experiments. 

keywords: Electromagnetic Corrections, Kaon Semileptonic Decay, Form 
Factors, Chiral Perturbation Theory. 
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1 Introduction 

Nowadays, it is widely recognized that vrvr scattering is the purest process al- 
lowing a direct probe of the Quantum Chromodynamics (QCD) ground state. 
This is due in part to the fact that pions are the lightest pseudo Goldstone 
bosons coupled to the vacuum thanks to the mechanism of spontaneous and 
explicit breaking of Chiral symmetry. Let Nf be the number of light quark 
flavors and denote by Fq the coupling in question, then, the Goldstone's, 0"^, 
satisfy, 

mi{m\v)) = i^'^'F.p,, a,6 = l,...,iV|-l, (1) 
where Af, are the usual axial currents. 
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The coupling Fq is given by the pion decay constant in the hmit of van- 
ishing masses for the Nj hght quarks, 



= Fo [1 + OK)] , q = u,d,s. (2) 

The non vanishing of Fq is a necessary and sufficient condition for the spon- 
taneous breaking of Chiral symmetry p. While the occurrence of the latter 
mechanism is a fact, the way it occurs remains uncertain. It is currently as- 
sumed that, resting on large- A^^c power counting grounds, the present mech- 
anism is due to important quark condensation in the vacuum. To confirm or 
reject this assumption one has to measure the quark condensate, {qq). To 
this end, it is convenient to consider the ratio [2], 

^ ^ (m„ + mrf)(gg) 
FfMl—' 

which vanishes for a vanishing quark condensate and equals one for large 
value of (qq). In order to fix the value for X from experiment, it is inter- 
esting to perform accurate measurements of observables that are sensitive 
to variations of X. This is the case of low-energy vrvr scattering. In fact, 
the latter is solely described in the threshold region in terms of the S-wave 
scattering lengths, 

1L_ [5a + 16/5 + 0{m,)] , (4) 



where, 



a = A-3X + 0{mq) , /3 = 1 + 0{mg) . (6) 

For instance, the isoscalar scattering length, Oq, varies by a factor 2 when X 
varies from 1 to zero. This sensitivity of vrvr scattering to the size of the quark 
condensate is not only a feature of the leading order but, on the contrary, it 
persists at higher orders allowing vrvr scattering to be the golden process for 
testing the mechanism of quark condensation f3]. 

The TTTT scattering is experimentally accessible in K^^ and pionium decays. 
In fact, the partial wave expansion of form factors displays the vrvr phase 
shifts as stated by the Watson final state interaction theorem jH The 
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data on phase shifts can then be translated into a model-independent deter- 
mination of scattering lengths by means of Roy equations [U]. Concerning 
the pionium characteristics such as its lifetime, r, and 2S — 2P strong energy 
level shift, AEs, they give direct access to im scattering lengths via |7|, 

r oc (a° - alf , AE, oc 2a° + al . (7) 

Once the results from the presently running DIRAC experiment are avail- 
able, Oq and Og should be determined with 5 to 10% accuracy 8 . On the 
other hand, the charged K^i decay has been measured by the E865 experi- 
ment [HI and the outgoing data have been analyzed in ^T] and ^2] 5 in- 
dependently. Before recalling the conclusions of both references, let us stress 
that the obtained value for Og is 7% accuracy and is compatible with the 
prediction of the standard power counting which rests on large Nc grounds. 
The analysis of Ref. ^T] relies on chiral symmetry constraints which corre- 
late the two scattering lengths. It yields a value for Oq consistent with the 
standard counting prediction. Furthermore, if one combines the extracted 
values for the scattering lengths, the ratio Q reads then, X ~ 0.94. As 
for the analysis performed in Ref. ^21; it combines the data of the E865 
experiment with an existing one in the isospin-two channel below 800 MeV 
without using the aforementioned chiral constraints. The conclusions of this 
analysis point out a discrepancy at the 1 — a level between the extracted 
value for Oq and the one predicted by the standard power counting. More- 
over, the corresponding value for the ratio © was found to be, X ~ 0.81. 
In view of the "disagreement" between the results of the two analysis, and 
before drawing off any conclusion about the size of the quark condensate, it 
is necessary to devote much more effort at both the experimental and the 
theoretical levels. In this direction, new precise measurements of charged and 
neutral Ke^ decay are currently taking data at CERN and FNAL 
respectively. This should be accompanied by an improvement of the accuracy 
in the theoretical prediction for the scattering lengths by evaluating isospin 
breaking effects in tttt scattering as well as in K^^ decays. While such effects 
in the former case are now under control at leading and next-to-leading or- 
ders [13 Cni UZl UHl Uni I2n|; we are interested in evaluating these effects at 
the same orders for the latter case. Recently, one of us has published the cal- 
culation of Isospin breaking effects in the K^^ decay of the neutral kaon [21] . 
The aim of the present work is to evaluate the same effects in the Ki^ decay 
of the charged kaon. 
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2 K^4: decay of the charged kaon 

The semileptonic Ki4, decay of the charged kaon is given schematically by, 

K+{p) ^ 7r+(pi) 7r-(p2) t{pi) u,{p,) , (8) 

where the lepton, i, is either a muon, /i, or an electron, e, and u stands for 
the corresponding neutrino. 

2.1 Matrix element 

The decay (jHl) is described in terms of an invariant decay amplitude, A'^~ , 
defined via the matrix element, 

= I (27r)"5(^) (pi +p^+pi+p^-p) [-1 A^-) , (9) 
with the on-shell conditions, 

p2 = Mi±, = MI^^pI = Ml^^p] = m^,pl = 0. (10) 

2.2 Form factors 

Let us introduce the notations, 

P = p^+p^^ Q = p^-p^^ L = pi+p^, N = pi-p^. (11) 

By Lorentz covariance, the decay amplitude can be parameterized in terms 
of three vectorial form factors, /"*"", g~^~ , r~^~, one anomalous form factor, 
h'^~ , and one tensorial form factor, T, 



GfV* 



Mk± 



In the preceding formula, a'^" is nothing else than the usual gamma matrix 
commutator. 



Vus denotes the Cabibbo-Kobayashi-Maskawa flavor- mixing matrix element, 
and Gp is the so-called Fermi coupling constant. Note that the normalization 
factors are written in powers of the charged kaon mass by pure convention. 
Let F~^~, G~^~ , -R"*"", and H~^~ , denote respectively the form factors, /"*"", 
g+- ^ r"*"", and h~^~ in the Isospin limit. At Tree-level, these quantities are 
given by the Current Algebra formulae [221 1221 1211 12^1 , 

Q+- 
Mk± 

V2Fo' 

Mk± f s^ + t^- 

2V2F0 V ; ' 

0. 

Herein, the Lorentz invariants, and si, represent the dipion and dilepton 
masses, respectively, 

S.= {P1+P2?, si = {pi+p,f. (12) 

The scalars, t.^ and Ut^, denote exchange energies between and the two 
pions, TT"*" and vr", respectively, 

U = (P-Plf, Mtt = {P-P2f ■ (13) 

2.3 Decay rate 

The decay rate F for process (jH)) is obtained from the differential decay rate, 

dV = —^d^y\A+-\\ (14) 
2Mk± ^11' ^ ^ 

spins 

by integrating over the differential phase space, 

= {271^6^'^^ {pi +P2+Pi+Pu-p) y< 

d^Pi <i^P2 d^Pi <i^Pu 
(27r)32Ei (27r)32E2 (27r)32Ez (27r)32|pJ ' 

with particle energies, 

= \IP1 + Ph (15) 



F+- 
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To this end, we will follow the approach used in [22] which consists on looking 
at decays as being two-body decays into a dipion of mass and a 
dilepton of mass si. The two systems subsequently decay in their own center- 
of-mass frames. To describe the decay distribution it is convenient to use, 
besides the invariant masses, s-^ and s/, the angles 9^,, 9i, and as illustrated 
in figure d Then, one has to express all of the scalar products obtained 
from momenta, pi, p2, Pi, and p^, or equivalently from, P, Q, L, and N, 
in terms of the five independent variables, s^, s/, 9i, and 0. Naturally, 
Isospin breaking effects enter the expression of these scalar products through 
the difference between the two pion masses. Since the mesonic final state in 
process (jHl) consists on two charged pions, one concludes that Isospin breaking 
does not affect scalar products, 

= s^, (16) 

Q2 = 4M2±-s,, (17) 

= si, (18) 

= 2mf-si, (19) 

P-Q = 0, (20) 

P-L = 1 (M|± - - s;) , (21) 

P-N = ^{Mli - - si) zi + {1 -zi)X cos 9i, (22) 

Q ■ L = Xcttt cos 9^^ , (23) 

Q-N = ziXa^ cos + ^ (1 - zi)a^ [ {Mj^± - - si) x 

cos 9t, cos 9i — sin 9^, sin 9i cos (p ] , (24) 

L-N = mf , (25) 
e^.p^L^'N'^pPQ'' = zi)Xa^y^ism9^sm9ism(P. (26) 

Herein, 

2 

zi = ^, 27 

Sl 

X = l\'/\Ml^,s^,si), (28) 



a, = Wl ^, (29) 
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and the function, 

A(x, z) = + if + — 2xy — 2xz — 2yz , (30) 

is the usual Kallen function. 

The same holds for the differential phase space which reads in terms of 
the five independent variables, 

d^^ = M^±N{s^, si)ds^dsid{cos6^)d{cosei)d(j) , (31) 

with, 

Nis^,s,) = (32) 

Obviously, the phase space is reached by integrating (p?T|l over the following 
range of variables, 

< < 27r 

< 9i < TT 

< < TT (33) 

mf < si < (Mk±-v^)' 

The last step consists on squaring the amplitude and summing over spins 
to get, 

^ 2Gl\V^s? 3,{s^,sue^,ei,<P) , (34) 

spins " 

with the following expression for the intensity spectrum, 

,2\ 



j5 = \f\'[{P-Lf-{P-Nf~s^{si~m] 
+ \g\^[{Q.Lf-{Q.Nf~Q\si- 



I I |2 2/ 2\ 

+ \r\ rrii [Si — rrii ) 

+ (r^7 + /^?*) \{P ■L){Q-L)-{P- N){Q ■ N)] 

+ (/*r + fr*)m1 [(P ■ L) - (P ■ N)] 

+ {rh + fh*)-^[{Q.N){P.Lf 

- {Q-L){P- L){P ■ N) - s^siiQ ■ N) + m]s^{Q ■ L 



+ {g*r + gr*)mf[{Q-L)-{Q-N)] 
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+ i9*h + gh*)j^[{P.L){Q-L){Q-N) 



- {P ■ N){Q ■ Lf + siQ\P ■ N) - m]Q\P ■ L 

- e.^.^L^N'^PPQ'^ [Ml4rg - fg*) 



Ml. 

- (P ■ N){rh - fh*) - (g ■ N){g*h - gh*) - m}{r*h - rh*)] 
+ {fT* + Tn,^^[s^{Q.N)-s^{Q-L)] 

+ {gT* + Tg*),^^[Q\P-L)-Q\P-N 

+ {rT* + Tr*) [{P ■ L){Q ■ N) ~ {P ■ N){Q ■ L)] 

+ {rT*-Tr*)^r^e,,p^L'^N''PPQ'^ 
2,Mk± 

- {hT* + Th*) [-Q\P . Lf + Q\P .L){P.N)- s^{Q ■ Lf 
+ s^iQ ■ L){Q ■ N) + s^siQ^ - s^Q^L ■ N 

- s^Q^N^ + 2Q\P-Ny + 2s^{Q-Nf]. (35) 

From the foregoing we conclude that, apart from the corrections affecting 
the form factors F~^~ , G~^~ , i?"*"" and H^^ directly, Isospin breaking en- 
ters the expression of the intensity spectrum through the additional terms 
proportional to the tensorial form factor, T. 



3 Form factors to one-loop level 

We present here a one-loop calculation of the decay amplitude for the 
process (jH)) including isospin breaking terms. The starting point is an effective 
chiral Lagrangian describing the dynamics of mesons, photons and leptons. 
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3.1 The effective Lagrangian 

In order to treat completely electromagnetic effects in Ki^ decays, not only 
the pseudoscalars but also the photon and the light leptons have to be 
included as dynamical degrees of freedom in an appropriate effective La- 
grangian [27]. The starting point is QCD in the limit m„ = = = 0. 
The resulting chiral symmetry, G = SU{3)l x SU{3)r, is spontaneously bro- 
ken to SU{3)v- The pseudoscalar mesons {n^K,!]) are nothing else than 
the corresponding Goldstone fields 0j = 1, . . . , 8) acting as coordinates of 
the coset space SU{3)l x SU{3)r/ SU{3)v- The transformation rules for the 
coset variables ul,r{4>) are 

gLULh{g,(j))~^ , 
UR{(j)) gRUnh{g,(j))-'^ , 
9 = {9l,9r) e Sf/(3)LxSf/(3)fl, (36) 

where h{g, 0) is the nonlinear realization of G j2Hl I2H]- 

As stated before, the photon field and the leptons ^, {i = e, /i) have 
to be dynamical. Thus, they most be introduced in the covariant derivative, 

= i [u^jiid^ - ir^)uR - u\{d^ - il^,)uL], (37) 

by adding appropriate terms to the usual external vector and axial-vector 
sources V^, A^. At the quark level, this procedure corresponds to the 
usual minimal coupling prescription in the case of electromagnetism, and 
to Cabibbo universality in the case of the charged weak currents, 

= v^ + a^-eQ'^^A^. (38) 

The 3x3 matrices Qf^R, Ql[ are spurion fields corresponding to electromag- 
netic and weak coupling, respectively. They transform as, 

QT^ - 9lQT^9{ , QT ^ 9RQT9i , (39) 

under the chiral group. In practical calculations, one identifies QJ^r with the 
quark charge matrix 

j 2/3 Q \ 
Q'^"^ = -1/3 , (40) 
\ -1/3 / 
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whereas the weak spurion is replaced by, 

Ql = -2^/2 Gi. ( | , (41) 


where Gp is the Fermi coupling constant and Vud, Vus are Kobayashi-Maskawa 
matrix elements. 

In order to take into account electromagnetic mass difference between 
charged and neutral mesons, it is convenient to define the following electro- 
magnetic and weak sources, 

QT'^ = u\QT-u,, QT = u^nQ^un (42) 

transforming as 

QT- ^ h{g,<P)QT'-h{g,<j>r\ 

QT ^ Hg^mTKg, <!>)-'. (43) 

With these building blocks, the lowest order effective Lagrangian takes 
the form 

- ^F^^F^''' + J2[I{i^ + e4-mi)£ + VeLiML], (44) 
e 

where ( ) denotes the trace in three-dimensional flavor space. 

The fow-energy constant Fq appearing in the preceding formula is an order 
parameter for chiral symmetry since it testifies to its spontaneous breaking. 
It represents the pion decay constant in the chiral limit, = rud = = 0, 
and in the absence of electroweak interactions. Explicit chiral symmetry 
breaking due to quark masses is included in, 

X+ = u\xul + u\x^ur . 
In practice, one makes the following substitution, 

ruu 

X ^ 2SoMq,ark = 2^0 | TUd | , (45) 

rUs 
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where Bq is an order parameter for chiral symmetry. It is related to the 
quark condensate in the chiral limit by, 

{qq) = ~F^Bo. (46) 

The low-energy constant Zq expresses explicit chiral symmetry breaking by 
electromagnetism. It is given by the electromagnetic mass of the pion as we 
will see below. 

In the absence of electroweak interactions, that is, for = rriii , a = 
0, ChPT is Isospin- invariant. In order to study Isospin breaking effects in 
ChPT processes, the usual chiral expansion in powers of p and is no more 
sufficient. One must also expand matrix elements in powers of the isospin 
breaking parameters, rrid — mu and a. On the other hand, the best accuracy 
ever reached in strong interaction observable measurements does not exceed 
the 5% level. Thus, we consider that an expansion to orders 0{md — mu) and 
0{a) is highly adequate for our purposes. Moreover, chiral expansion and 
Isospin breaking expansion have to be related in a consistent way in order 
to obtain reliable results. We adopt an expansion scheme where the Isospin 
breaking parameters are considered as quantities of order in the chiral 
counting, 

0(m, - m„) = Oia) = OK) = 0(/) . (47) 

Therefore, tree level calculation corresponding to leading chiral order is char- 
acterized by chiral orders cited in (j47|) . Concerning one-loop level calculation 
which corresponds to next-to-leading chiral order, it is characterized by the 
following chiral orders. 



O {p\md - m„)) , O (m.K - ^«)) , 0{p''a) , 0{mga) . (48) 

We have now all necessary elements to calculate any Green function in 
the framework of ChPT including Isospin breaking effects. For example, the 
leading chiral order expressions for light meson masses are found to be, 

M^o = Ml, (49) 
Ml± = Ml + 2Zoe^F^, (50) 

Mlo = Ml + ^{Ml- Ml) , (51) 
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= Ml-^{Ml-M'^)+2Zoe'F^, (52) 

= 1 (4Mi - Ml) . (53) 

Herein, M^r and Mk represent respectively pion and kaon masses in the 
absence of isospin breaking, 

Ml = 2Bom, Ml = Boim + nis) , 2m = m„ + , (54) 

e measures the rate of SU{2) to SU{?>) breaking. 



4 m, — m 



(55) 



At next-to-leading chiral order, one-loop calculation is involved. Vertices are 
extracted from Lagrangian (jl^ . Meson masses in the propagators as well 
as in the vertices can be identified with expressions P9|) - (f53|l . As is well 
known, loops are ultraviolet divergent. To remove divergences, renormaliza- 
tion should be employed. The procedure consists on adding to Lagrangian 
(HD) suitable counter-terms [201 EH E2 IMl I2I| generating exactly the same 
divergences but with opposite sign. Moreover, the cancellation should occur 
order by order in the chiral expansion as dictated by renormalizability prin- 
ciples of effective field theories. Counter-terms are modulated by low-energy 
constants which are order parameters for chiral symmetry. In order to deter- 
mine these constants one proceeds as follows. Let C be either a low-energy 
constant or a combination of low-energy constants and T an observable very 
sensitive to variations of C. One first calculate the expression of T in the 
framework of ChPT to any given order and then match the obtained ex- 
pression with an experimental measurement of T. It is clear that the value 
for C deduced from this matching does not constitute a genuine determina- 
tion of the low-energy constant. In fact, it represents the value of C at the 
given chiral order and with the accuracy of the experimental measurement. 
Note that the method of effective Lagrangian has the disadvantage of an 
infinitely increasing number of low-energy constants when going to higher 
and higher orders in the low-energy expansion. For instance, two constants 
in the strong sector, Bq and Fq, and one constant in the electroweak meson 
sector, Zq, parameterize the leading chiral order. At next-to- leading chiral 
order, one has ten low-energy constants in the strong sector, Li , . . . , Lio, 
fourteen constants in the electroweak meson sector, Ki , . . . , K14, and seven 
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constants in the electroweak leptonic sector, Xi , . . . ,Xj. The constants, 
Lj, Ki and Xj, are divergent. They absorb the divergence of loops via the 
renor mahzat ion , 

U = L[(/i) + r,A, ^ = 1,... ,10, (56) 
= Kl{fi) + J:,X, I = 1,... ,14, (57) 
X, = X[(/i) + S,A, I = 1,... ,7. (58) 

Herein, A corresponds to pole subtraction in the MS dimensional regulariza- 
tion scheme (see appendix). The beta-functions, Fj, Ej and Sj, can be found 
in jSO], IHI] and |2Z|, respectively. The scale /x cancels in observables as can 
be seen from the renormalization group equations, 

^K/^2) = ^K/^i) + Y^ln^, (59) 
and similar for Ki and Xj. 



3.2 Leading order 

Let us start with the tree-level calculation corresponding to the leading chiral 
order. The different topologies for K^/^ decays at tree level in perturbation 
theory are drawn in figure |21 Decay amplitudes for the different processes 
were calculated in where the expressions for the corrected form factors 
in the charged channel were found to be, 

r- = F^~, 
9^- = 



R+- + AZoe'F^ 



2 r2 



y^Fo Si - M|,± 



2e^F^ [ t^-u^ Q-L + Q-N \ Mk± 



T 



si - Ml^ M|± - -2p-pJ V2Fo ' 
Ae'^F^ miMK± Mk± 



3.3 Born contribution 

We proceed with the calculation of form factors and g~^~ to next-to- 
leading chiral order. Feynman diagrams representing the amplitude will be 
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separated into two sets: photonic and non photonic diagrams. The non pho- 
tonic set is drawn in figure El The calculation of these diagrams is standard 
in field theory. The starting point is Lagrangian ()44j) . For the non linear 
realization of chiral symmetry we will use the exponential parametrization, 

ur = u{ = expj^l , (60) 

where $ is the linear realization of SU{'i) and can be decomposed in the 
basis of Gellmann-Low matrices as, 



a=l 



In terms of physical fields the matrix $ can be written, 



$12 = , $13 = -V2K+ , $21 = V2- 

$22 = (-l + ^=\'n'+U+^\v. 



$23 = -V2K\ $31 = V2K-, $32 = -V2K^, 

$33 = -^{i2n' + v) . (61) 

The two mixing angles ei and €2 relate 0^ and 0^ to the mass eigenstates 7r° 
and 7], 

vrO \ / 1 ei ^ / 

V J V -e~2 1 

Notice that ei = £2 = e at leading order. The next step consists on expand- 
ing Lagrangian to fifth order in pseudoscalar fields, generating Feyn- 
man rules and drawing allowed topologies. We then calculate pseudoscalar 
propagators and derive masses and wave function renormalization constants. 
Finally, we expand the next-to-leading order Lagrangian to third order in 
pseudoscalar fields and obtain the counterterm contribution. 

Let us denote by 6F and 6G the next-to-leading order corrections to the 
F^^ and form factors, respectively. 



(62) 
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+- 



M 



V2Fo 



1+6G 



The expressions for 6F and 6G are lengthy. Therefore we will separate them 
to different contributions depending on the topology of the Feynman diagram 
representing a given contribution. Notice that in the following we will use 
the same notations as in reference ^21j| where definitions and expressions for 
loop integrals have been given in the appendix. 

Let us start with Born contribution. This contribution is obtained from 
diagram (a) in figure El We take the corresponding vertex from Lagrangian 
()44|) and multiply by the wave function renormalization constant factor to 
obtain, 



6F 



+ 



6G 



3(3 + -^ ) Ao(M^o; 



3 1 



+10Ao(M^±) + 6Ao{Mko) + 8Ao{Mk±)] 



2 



4 



1 



Ml ' Ml ^' m2 



Aq{M^) 2^oK) 

I 



— ^ 9 + 21n-i- + 41n-| + 21n^ 
167r^ \ M^ M^ mf 



6 



{2AKi + 247^2 + 207^5 + 207^6 + SXg 



(63) 



The Born contribution is infrared divergent. This divergence emerges from 
the wave function renormalization constants of charged kaon, pions, and 
lepton. 

3.4 Counter-terms contribution 

The contribution in question follows from diagram (a) in figure |21 with the 
help of the next-to-leading order Lagrangian and reads, 



(5F = -2 [32pi ■ P2L1 + 4(M|.± + 2Ml^ - + 2pi ■ - 2pi ■ Pu)L2 
^0 

+ (M|.± + 2M2± - + lOpi ■p2 + 2p-pi-2p-p2- 2pi ■ p,)L^ 
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+ 4(2M| + hMl)U + 2(M|± - + 3M2o)L5 
+ (mf + ■pi.)L9] 

+ — {l2Ki + 847^2 + 19/^5 + 377^6 + 97^12 - 30Xi) , 
9 
2 

^0 

- (M|± + 2M2± - + 2pi ■p2 + 2p-pi-2p-p2- 2pi ■ p^)L^ 

+ A{Ml + 2M|)L4 + 2(M|± - M^^ + 3M2o)L5 + (m^ + 2p; ■ p,)L^] 

+ ^(127^1 + 127^2 

+36^3 + 187^4 + 7i^5 + 257^6 + ^Ki2 + 6Xi) . (64) 



3.5 Tadpole contribution 

Tadpoles are shown in diagram (c) of figure |21 and contribute to the form 
factors by the following, 

= T7^[^o(M^o) + 9Ao(M,) 
12-^0 

+8Ao(M^o) + 4Ao(M,±) + 8Ao(M^±)] , 
5G = ^ [Ao(M,o) + Ao(M,) + 4Ao(M^±) + 4Ao(M^±)] . (65) 

3.6 The s-channel contribution 

The remaining contribution to the non photonic part of the decay ampli- 
tude comes from loop diagrams with two pseudoscalar propagators. This 
two-point function contribution will be separated in three parts depending 
on the Lorentz scalar governing its underlying kinematics. The s-channel 
contribution comes from diagram (c) in figure 121 and reads, 

6e 

+ _ Ao(M^) - 4Ao(M.±) - 2Ao(M^o) - AA^{Mk±) 
+ 3 [2M2± - M^o + 2pi • p2 
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6G 



+-|(M2 + 2pi-p2) 



+ 31 



2e 

7! 



M^oBo{-pi-p2,Mr„Mr, 



Vie 



- -j={M^ + 3pi-p2)Bo{-pi-p2,M^,Mr,) 

- 8(4M2± - 3M^o + pi ■ P2)Bi{-pi - p2, M^±,M^±) 

- 4(M2± + p, . p2)B,{-p, - P2, Mko, Mko) 

- 8(4M2± - 3Mlo + pi ■ P2)Bi{-pi - p2, Mk±, Mk±) 
+ 45oo(-pi -P2,M,±,M,±) 

+ 25oo(-pi -p2,Mxo,M^o) 

+ 45oo(-pi -p2,Mx±,M^±) 

+ 8(M2± + . p^)Bu{-pi - P2, M^±, M^±) 

+ 4(M^± +pi-P2)Bu{-Pi-P2,Mko,Mko) 

+ 8(M2± +pi ■p2)5ii(-pi -P2,M;^±,M^±)} , 

= -l^{2Boo{-Pi-P2,M^±,M^±) 

^0 

-5oo(-pi -p2,M^o,M^o) + 25oo(-pi -P2,Mk±,Mk±)} 



3.7 The ^-channel contribution 

Diagram (d) in figure El generates the somewhat lengthy t-channel contribu- 
tion, 



SF 



1 



12F2 



{-6AoiMK<^) - 2AoiMK±] 



+ Q[Ml^-Mlo-p-pi 



2e 



+ ^[Ml-Ml-2p-p{) 



Boipi-p,M^o,MKo) 

- 3 [3M^o - Ml - 2Ml± +2p-pi 

--^ {Ml - M^ + 8p ■ pi) 5o(pi - p, M^, M;^o) 

- 6 [M|.± + - M^o - 2p ■ pi 
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B,{p,~p,M^o,Mko) 



+ 3 [2Ml± - 2M2± - 3M^o + Ml 



+^(M^ + Ml-p.p,) 



Bi{pi -p, Mr,, Mko) 



- 4(4M2± - 3M^o - p ■ Pi)B,{p, - p, M^±,Mk±) 

+ 12 (^1 - ^ Booipi - P, M^o, Mko) 

+ 12(2+^^ Boo{pi-p,Mr„MKo) 

+ 20Booipi -P,M^±,Mk±) 



--^(2Mi-M,2-p.p,) 
+ 6[M|±-2M2±+p.pi 
+ -^(2Mi-M2-p.pO 



Bu{pi-P,M^o,Mko) 



Bnipi - P, M^, Mko) 



+ 4(2M|.±-M2±-p-pi)Sn(pi-p,M^±,Mx±)} , 
SG = -^{-6Ao{Mko)-2Ao{Mk±) 



+ 6[M^i-M^o-p-pi 



2e 



+^{Ml-Ml-2p-p,) 



Bo{pi-P,M^o,Mko) 
3 [3M^o - Ml - 2Ml± +2p-pi 
--^{M^-Ml + 8p-p,) Bo{pi-p,Mr„MKo) 
6 [M|^± + M^± - M^o - 2p • pi 

5i(pi-p,M^o,M;^o) 



+ ^(2M^-P-Pi) 



+ 3 [2M|-± - 2M^± - 3M^o + M, 



4e 



+ ^iM^ + Ml-p-pi) 



Bi{pi-p,Mr„MKo) 
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- 4(4M2± - SM^o -p.p,)B,{p,- p, M^± , Mk± ) 
12e 

+ -^Boo{pi-P,M^,Mk) 



12 ( 1 + -^J Soo(pi - p, M^, M;^o) 
- 4Boo{pi-P,M^±,Mk±) 



-^(2Afi-M,^-p-pO 
+ 6 [M|.± - 2M2± + p ■ pi 
+ ^(2Mi-M^2-p.pO 



Bu{pi-P,M^o,Mko) 



5n(pi-p,M^,Mxo) 
+ 4(2M|:± - - p ■ pi)Bu{pi - P, M,±, M^±)} . 



3.8 The w-channel contribution 

Finally, the w-channel contribution follows from diagram (e) in figure 01 



6F 



6G 
1 

6^ 



{Ao(M,±)+Ao(Mx±; 



+ 6(M2± - M^o + p ■ P2)Bo{p2 - p, M^±,Mk±)} , 
3.9 Soft virtual photon contribution 

The various topologies of Feynman diagrams containing a virtual photon are 
drawn in figures 111013 and[7| The contribution of each figure will be given 
separately. 

Figure SI contribution 
diagram 4. (a) 



6F = -^{2Bo{-p,0,Mk) + B,{-p,0,Mk)} , 6G 



diagram 4. (b) 



SF = 6G = 0. 



(66) 
(67) 
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diagram 4. (c) 

SF = 5G = {2Boi-pu 0, M^) + B^i-p,, 0, M^)} . (68) 

diagram 4. (d) 

SF = -5G =^ {25o(-p2, 0, M^) + 5i(-p2, 0, M^)} . (69) 

Figure [H contribution 
diagram 5. (a) 

(5F = 6G 

= -e^[Ap-p2Goi-p2,-p,fn^,M^,MK) 

+ 5o(-p2, 0, M^) + 5o(-p, 0, Mk) - Boip2 - P, M^, Mk)] . (70) 

diagram 5. (b) 

= [4p-piCo(-pi,-p,m^,M,,Mx) 

+ Bo{-pi, 0, M^) + Boi-p, 0, Mk) - 5o(pi - p, M,, M^) 

+ 4p ■ PiCi(-pi, -p, 0, M,, Mk) + 8p ■piC2(-pi, -p, 0, M,, M^) 

+ fii(-pi,0,MO + 25i(-p,0,MK) 

+ Bo{pi-p,M^,Mk)-Bi{pi-p,M^,Mk)] 

5G = [4p-piCo(-pi,-p,m^,M,,MK) 

+ Bo{-pi, 0, M,) + Boi-p, 0, Mk) - 5o(pi - p, M^, Mk) 

+ Bo{pi-p,M^,Mk) + B^{p^-p,M^,Mk) 

+ fii(-pi, 0, M^ + Ap ■ piGii-pu -p, 0, M,, Mk)] . (71) 

diagram 5. (c) 

5F = ^[6Bo{-pi,0,mi) 

+ Bo{-p,0,Mk) + 2B,{-p,0,Mk) 

+ 12p ■ piGo{-pi, -p, m^, mi, Mk) + 2mfGi{-pi, -p, 0, mi, Mk) 
- %{Ml - 2p ■ pi)G2i-pi, -p, 0, mi, Mk)] , 
e2 

6G = -[35o(-p,0,Mk) 
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+ 12p- piCo{-pi,-p,'mj,mi,MK) 

- 6{inf -2p-pi)Ci{-pi,-p,0,mi,MK) 

+ 6M]^C2{-pi,-p,0,mi,MK)] . 

diagram 5. (d) 

SF = -e2[-5o(pi,0,MO-5o(-p2,0,M^) 

- 5i(pi,0,M^) + 5i(-p2,0,M,)-2Bi(-pi-p2,M,,M,) 
+ 4pi ■ P2Co{pi, -p2, m^, M^, M^) 

+ 4pi ■p2Ci(pi,-p2,0,M,,M,) 

- 4pi ■ p2C2(pi, -P2, 0, M,, M,)] , 

(JG = -e2[-5o(pi,0,M.)-5o(-p2,0,M,) 

- Si(pi,0,M,) -Si(-p2,0,MO 
+ Api ■ P2Co{pi, -p2, m^, M^, M^) 
+ 4pi ■p2Ci(pi,-p2,0,M^,M^) 
+ 4pi ■ p2C2(pi, -P2, 0, M^, M^)] . 

diagram 5. (e) 

5F = SG 

= -^[-^Bo(pi,0,mi) 

- 5o(-pi,0,M,)-25i(-pi,0,M,) 
+ 12pi • PiCq{pu -Pu m^: rni: M^) 

- 2mfCi{pi,-pi,0,mi,M^) 

+ 6(M^ + 2pi • pOC2(Pi, -Pi, 0, m,, M^)] . 

diagram 5. (f) 

SF = -^[-6Bo{pi,0,mi) 

- Bo{-p2, 0, M^) + 65o(-p2 - P/, m,, M„) + Si(-p2, 0, M^) 
+ 12p2 • PiCoiPi, -P2, m^, mi, M^) 

- 2mfCi{pi,-p2,0,mi,M^) 

+ 6M^C2ipi,-p2,0,mi,M^)] , 

SG = -j[-6Bo{pi,Q,mi) 
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- 5Bo(-P2, 0, M,) + 65o(-p2 - Pi, mi, M^) - Si(-p2, 0, M,) 

+ 12^2 ■ PiCoipi, -P2, m^, mi, M^) 

+ 2m^Ci{pi,-p2,0,mi,M^) 

+ QM^C2ipi,-p2,0,mi,M^)] . (75) 

Figure |H1 contribution 
diagram 6. (a) 

5F = -^[-2Bo{-pi-p,,0,Mk)-2Bo{-p,0,Mk)- B^{-p,0,Mk) 
o 

+ 12(M^ +P1-P2-P-P1+P- P2)Co{-p, -pi - pu, 0, Mk, Mk) 
+ Q{Ml +pi-p2-3p-pi + 3p- P2)Ci{-p, -pi - p„ 0, Mk, Mk) 

- 12{p -pi-p- P2)C2i-p, -pi - Pu, 0, Mk, Mk) 

- Q{p-pi-p- p2)Cu{-p, -pi - pu, 0, Mk, Mk) 

- 6{p-pi-p- P2)Ci2{-p, -pi - Pu, 0, Mk, Mk)] , 

6G = e^Cooi-p,-pi-p„0,MK,MK). (76) 



diagram 6. (b) 



SF = 50 = Bo{-pi - Pu, 0, Mk) • (77) 



diagram 6. (c) 



6F = 6G 
4e2 



— [Boi~pi-p„0,MK) 

m]Ci{-pi - Pu, -pi, 0, Mk, mi) 

mfC2{-pi - Pv, -pi, 0, Mk, mi)] . (78) 



diagram 6. (d) 

6F = 6G 



— [ABoi-pi - Pu, 0, Mk) + 4So(pi, 0, M^) + 2Bi{pi, 0, M^) 

+ Q{Ml + 2Ml-m^i 

+2pi ■p2-2p-pi + 2p-p2- 2pi ■ Pu)Cq{pi, -pi - Pu, 0, M^, Mk] 
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+ 3(M| + - + 2pi ■ p2 

-2p ■pi + 2p-p2- 2pi ■ Pu)Ci{pi, -pi - p,, 0, M,, Mk)] . (79) 

diagram 6. (e) 

5F = ^[2Bo{p2,0,M^)-ABoi-pi-p,,0,MK) + B,ip2,0,M^) 



+ 3(M|. + 2M2-m2 

+2pi ■p2-2p-pi + 2p-p2- 2pi ■ Pu)Co{p2, -pi - Pv, 0, M^, Mk) 
+ 18p-p2Ci(p2,-p« -p.,0,M,,Mx) 

- &{Ml - 2Ml + - 2pi ■ p2 + 2pi ■ p,)C2ip2, -pi - Pu, 0, M^, M,^) 
+ 12CooiP2, ~Pi-Pu,0,M^,Mk) 

+ 6p ■ p2Cii{p2,-pi - Pu,0,M^,Mk) 

- ?,{Ml-2Ml + m] 

-2pi ■ p2 + 2pi ■ Pu)Ci2{p2, -pi - Pu, 0, M^, Mk)] , 

5G = --[2Boip2,0,M^)-ABoi-pi~Pu,0,MK) + Bi{p2,0,M^) 

+ 3(M|. + 2M2-mf 

+2pi ■ p2 - 2p ■ pi + 2p ■ p2 - 2pi ■ Pu)Co{p2, -pi - Pu, 0, M^, Mx) 
+ 18p-p2C^i(p2,-p« -p.,0,M,,Mj^) 

- 6(M|- - 2M2 + - 2pi ■ p2 + ■ p.)C2(p2, -Pi - P., 0, M^, M^^) 
+ 6p-P2Cu{p2,-Pi-Pu,0,M„,Mk) 

- ?,{Ml-2Ml + m] 

-2pi ■ p2 + 2pi ■ Pu)Ci2{p2, -pi - Pu, 0, M^, Mk)] ■ 

Figure [7| contribution 
diagram 7. (a) 

SF = ^{-2Bo{-pi-Pu,0,MK)-2Boi-p,0,MK)-B,{-p,0,MK) 
o 

+ 12(M2 + . p2 - p . Pi + p . p2)C^o(-p, -Pi - p., 0, M;^, Mk) 
- 12{p -pi-p- P2)Ci{-p, -pi - Pu, 0, Mk, Mk) 
+ Q{Ml +P1-P2-P-P1+P- P2)Ci{-p, -pi - Pu, 0, Mk, Mk) 
+ 2mfCi{-pi,-pi-pu,0,mi,MK) 
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- 12(p • pi - p • P2)C2(-p, -Pi - P., 0, Mk, Mk) 
+ 2mfC2{-pi, -Pi - Pu, 0, mu Mk) 

+ 2m'{C2{-p,-pi,0,MK,mi) 

- 6{p-pi-p- P2)Cn{-p, -pi - p,,, 0, Mk, Mk) 

- 6(p • pi - p • p2)Ci2(-p, -pz - Pz., 0, Mk, Mk) 

- 12mf(M^+pi •p2 -P-P1+P-P2) X 

D2{-p, -Pi, -Pi - Pu, 0, Mk, mi, Mk) 

- 12mf{M^+pi-p2-p-pi+p-p2) X 







-P: - 


-Ph -Pi -Pu,0,Mk, 


tuuMk) 


+ 


12mf (pi 


■Pl - 


- P2 • Pl) X 






^22 


i-P, 


-Pi, -Pi -Pu,0,Mk 


,mi,MK) 


+ 


12mf{p ■ 


Pl - 


• p ■ P2) X 






D33 


{-P, 


-Ph -Pi -Pu,0,Mk 


,mi,MK) 


+ 


12mf (p ■ 


Pl - 


•P-P2) X 






D12 


{-P, 


-Pi, -Pi -Pv,^,Mk 


,mi,MK) 


+ 


12mf{p ■ 


Pl - 


p ■ P2) X 








{-P, 


-Ph -Pi -Pu,0,Mk 


,mi,MK) 


+ 


12mf{p ■ 


Pl - 


■ P ■ P2 + Pl ■ PZ - P2 ■ 


Pl) X 




D23 


{-P, 


-Ph -Pi -Pu,0,Mk 


,mi,MK)} , 




-e' [Coo 


{-P, 


-Pl - p., 0, Mk, Mk 


) 


+ 


2m^Doo{ 




-Pi, -Pi -Pu,0,Mk, 


mi,MK)] . 



diagram 7. (b) 

5F = 5G 

= f^{25o(-pi-p,,0,MK) 

+ 2fio(pi,0,MO + Bi(pi,0,MO 

+ 12p-P2Co{pu-pi-p„0,M^,Mk) 

+ 6p-P2Ci{pi,-pi-Pu,0,M^,Mk) 

- 2mfCi{-pi,-pi - p^,0,mi,MK) 

- 2mfC2i-pi,-pi - p^,0,mi,MK) 

- 2mfC2{pi,-pi,0,M^,mi) 
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- 12mfp ■ P2-D2(pi, -Pi, -Pi - Pu, 0, M^, mi, Mk) 

- 12mfp • P2-D3(pi, -Ph -Pi - Pu, 0, M^, mu Mk)] ■ 
diagram 7. (c) 

5F = -^{-4Bo(P2,0,M,) 

+ 2So(-pi-p.,0,MK)-2Si(p2,0,M^) 

- mi-mt 

-2p ■ p2 - 2pi ■ Pu)Co{p2, -pi - Pu, 0, M^, Mk) 

- 3(Mi-m^ 

-Api ■p2-2p-p2- 2pi ■ Pu)Ci{p2, -pi - Pu, 0, Mtt, Mk 

- 2m'{Ci{-pi,-pi- p^,0,mi,MK) 

+ 12{MI +P1-P2-P- Pi)C2{p2, -Pi - P., 0, M^, Mk) 

+ AmfC2ip2,-pi,0,M^,mi) 

- 2mfC2i-pi,-pi-p,„0,mi,MK) 
+ 6CooiP2,-pi-p.,0,M^,MK) 

+ Gpi-p2Cn{p2,-pi-Pu,0,M^,MK) 

+ &{Ml +P1-P2-P- Pl)Ci2ip2, -Pi - Pu: 0, M^, Mk) 

+ 6mf{MK -mf -2p-p2-2pi-p^) X 

D2{P2: -Pi, -Pi - Pu, 0, M^, mi, Mk) 
+ 6mf{Ml -mf -2p-p2-2pi-p^) X 
D3{P2, -Pi, -Pi - Pu, 0, M^, mi, Mk) 
+ 12mfpi ■ PiD22iP2, -Pi, -Pi -Pu,0,M^,mi,MK) 

- l2mf{Ml+pi-p2-p-pi)x 

D33{P2, -Pi, -Pi - Pu, 0, M^, mi, Mk) 

- 12mfpi ■ P2Di2{p2, -Pi, -Pi - Pu, 0, M^, mi, Mk) 

- 12m]p, . p2DniP2, -Pi, -pi Pu, 0, M^, mi, Mk) 

- 12m^{M^ + pi ■ p2 - p ■ pi - pi ■ pi) X 

D23{P2, -Pi, -Pi - Pu, 0, M^, mi, Mk)} , 

SG = -{-ABo{p2,0,M^) 

+ 2Bo{-pi-p„0,MK)-2Biip2,0,M^) 

- mi-m^i 
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-2p • p2 - 2pi ■ p^)Co(p2, -Pi - Pu, 0, M^, Mk) 

- ml -ml 

-4pi • p2 - 2p • p2 - 2pi • Pv)Ci{p2, -pi - Pu, 0, M^, Mk) 

- 2rnfCi{-pi,-pi-p^,0,mi,MK) 

+ l2{Ml +PI-P2-P- Pi)C2{p2, -Pi - Pu, 0, M^, Mk) 

+ AmfC2{p2:-PhO,M^,mi) 

- 2mfC2{-pi,-pi-Pu,0,mi,MK) 

- 6Coo(p2, -Pi - Pu, 0, Mtt, Mk) 

+ 6pi ■ P2Cii{p2, -pi - Pu, 0, M^, Mk) 

+ Q{Ml +P1-P2-P- Pi)Ci2(p2, -Pi - Pu, 0, M^, Mx) 

+ 6mf{M]^ -w?i -2p-p2- 2pi ■ p^) x 

^2(P2, -Pi, -Pi - Pv, 0, M^, m,, Mk) 
+ 6mf{Ml -mf -2p-p2-2pi-p^) X 

D^{p2, -pu -Pi - Pu, 0, M^, mi, Mk) 
+ 12mfpi ■ piD22{p2: -Pi, -Pi - Pu, 0, M^, mi, Mk) 

- 12mf{M^+pi-p2-p-pi) X 

D33{P2, -Ph -Pi - Pu, 0, M^, mi, Mk) 

- 12mfpi ■ P2Di2{p2, -Pi, -Pi - Pu, 0, M^, mi, Mk) 

- Um^pi ■ p2Di3{p2, -pi, -pi - Pu, 0, M^, mi, Mk) 

- l2m]{Ml+pi-p2~p-pi-pi-pi)x 

D23{P2, -Pi, -Pi - Pu, 0, M,, mi, Mk)} ■ (82) 

4 Soft photon bremsstrahlung 

Virtual photon corrections to Kj(4 decay rate generate infrared divergencies. 
These cancel, order by order in perturbation theory, with the ones coming 
from real bremsstrahlung corrections. Assume that the emitted photons are 
soft, that is, their energies are smaller than any detector resolution, uj. It 
follows that radiative and non-radiative decays cannot be distinguished ex- 
perimentally and emission of real soft photons should be taken into account. 
Note however that only single soft photon radiation is needed to one-loop 
accuracy. 
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4.1 The decay amplitude 



A general feature of photon bremsstrahlung is that, in the soft photon ap- 
proximation, the bremsstrahlung amplitude is proportional to the Born am- 
plitude. Since we deal only with isospin breaking corrections to the F and 
G form factors, the Born amplitude is taken, all along this section, to be, 

= i^GFv:MPu)i,{i-i>{pi){P' + Qn- 



The contribution of form factors F and G to the Bremsstrahlung amplitude 
can be read off from diagrams in figure |H1 

Let e and q be, respectively, the polarization vector and the momentum of 
the radiated photon. The evaluation of diagrams in figurelHlis straightforward 
and read, to first order in the photon energy. 



A 



+ -7 



eAi 



p-e* , pi-e 



p-q 



+ 



+ 



P2 



Squaring the matrix element 
obtain 



(83) 

pi-q pi-q P2-q • 
and summing over polarizations, we 



pol. 



-e^\Ai 



Ml 



X 



{p ■ qf 
2p-pi 



+ 



mf 



+ 



Ml 



+ 



{p-q){pi -q) 
2pi ■ Pi 



{pi ■ (pi ■ 

{p-q){p2-q) 
2p2 ■ Pi 



+ 



Ml 



+ 



(P2 ■ qf 
'^P-Pi 
{p-q){pi ■ q) 
2pi ■ P2 



{Pi ■ q) {Pi ■ q) {P2 ■ q) (pi ■ q) {Pi ■ q) {P2 ■ q) . 

The preceding expression is singular for vanishing momentum of the soft 
photon. We shall attribute a small but non-vanishing mass to the photon, 
m^, in order to regularize this singularity. 



4.2 The decay rate 

The Ki4^^ differential decay rate is obtained by squaring the matrix element 
(jSni), summing over spins and polarizations and integrating over the photon 
phase space, 

d^^ = (27r)'^5(^) {pi+P2+Pi + Pu + q-p) X 
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d^Pi d^P2 d^Pi d^Pv 



(27r)32Ei (27r)32E2 {2tt)'^2Ei {2n)^2\pJ (27r)32|q| ' 

Using the definition of bremsstrafilung integrals as given in tlie appendix of 
reference the Ki^^ differential decay rate takes the following form in the 
soft photon approximation, 



dT, 



2M 



K 



+-|2 
B I 



spins 



+ I{pi,pi,m^,uj) + I{pi,pi,m^,uj) + I{p2,P2,m^,uj) 
- 21{p,pi,m^,uj) + 21{p,p2,m^,uj) - 21{p,pi,m^,uj) 

+ 2/(pi, -pi, m^, u) - 21{p2, -pi, m^, u) - 2J(pi, -p2, m^, uj) 
where non-singular terms have been dropped out. 



4.3 Cancellation of infrared divergencies 

The infrared divergent part of Ki^^ differential decay rate can be extracted 
from (jH^ using the definition of bremsstrahlung integrals from the appendix 
of reference 



dT 



IR 



spins 



X 



2 + pi-p2T{pi,-p2,M^,M^) 

-P ■ Pit{p, Pi, Mk, M^) + p ■ P2r{p, p2, Mk, M^) 
-p ■ piT{p,pi,MK,mi) -pi ■ piT{pi,-pi,M^,mi) 



+P2- piT{p2,-pi,M^,mi) 



In 



On the other hand, the infrared divergence coming from virtual photon cor- 
rections to F and G form factors only can be read off from table ^ and is 
denoted by dV^^. It is easy then to check that infrared divergencies cancel 
at the level of differential decay rates. 



dT^^ + dT 



IR 



0. 



^5) 
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5 Perspectives 



In this work we studied the decay process, — > ■n^Ti'l^v^, taking into 
account Isospin breaking effects. These come mainly from electroweak inter- 
actions and generate corrections proportional to the fine structure constant, 
a, and to the difference between up and down quark masses, m„ — m^. 

The interest in this decay comes from the fact that the partial wave ex- 
pansion of the corresponding form factors involves tttt scattering phase shifts. 
The latter can be related in a model-independent way to the tttt scattering 
lengths which are sensitive to the value of the quark condensate. Thus, a 
precise measurement of form factors should allow accurate determination of 
scattering lengths and, consequently, give precious information about the 
QCD vacuum structure. 

Scattering lengths are strong interaction quantities. On the other hand, 
any K^^^ decay measurement contains contributions from all possible inter- 
actions, in particular, from electroweak ones. Therefore, it is primordial to 
have under control Isospin breaking effects in order to disentangle the strong 
interaction contribution from the measured form factors. The present work 
was guided by this motivation and, to this end, analytic expressions for form 
factors were obtained including Isospin breaking effects. These expressions 
are ultraviolet finite, scale independent, but infrared divergent. We showed 
that this divergence cancels out at the differential decay rate level if we take 
into account real soft photon emission. 

Our work should be completed by, 

• a parametrization of form factors in the presence of Isospin breaking, 

• a full treatment of the radiative decay, K^^^. 
Calculations in this direction are in progress. 
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Figure 1: Angles and unit vectors used in the kinematical description of 
decays. E^ttt and E;^,; are the planes defined in the kaon rest frame by 
the pion pair and the lepton pair, respectively. {9i), the angle formed 
by Pi (pi) , in the dipion (dilepton) rest frame, and the line of flight of the 
dipion (dilepton) as defined in the kaon rest frame. (/>, the angle between the 
normals to Stttt and T^i^^. v is a unit vector along the direction of flight of 
the dipion in the kaon rest frame, c (d) is a unit vector along the projection 
of Pi {pi) perpendicular to i). 
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(c) (d) (e) 

Figure 2: Feynman diagrams representing K^^ decay amplitudes at tree level. 
Wavy lines stand for photons. Only diagrams (a) and (b) contribute to the 
decay amplitudes J\P^ and J\^~ corresponding to K(,i decays of the charged 
kaon to neutral pions and of the neutral kaon to a neutral and a charged 
pion, respectively. 
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(d) 



Figure 3: Feynman diagrams representing the K£4^ decay amplitude of tlie 
charged kaon at one-loop. Is shown only the non photonic topology. Diagram 
(a) represents Born and counter-terms contributions. Tadpole contribution 
is given by diagram (b). Diagrams (c), (d) and (e) stand for contributions 
from the s-, t- and ti-channels, respectively. 
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(b) (c) (d) 



Figure 4: Virtual photons in Feynman diagrams for the K^^ decay in the 
charged channel. The topology consists on attaching a photon on the pure 
vertex from one side and a leg from the other side. 
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(d) (e) (f) 

Figure 5: Virtual photons in Feynman diagrams for the K£4 decay in the 
charged channel. The topology consists on attaching a photon between two 
legs of the pure vertex. 
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(d) 



Figure 6: Virtual photons in Feynman diagrams for the K(^^ decay in the 
charged channel. One takes the pure vertex mediated by a K'^ pole and then 
attaches a photon between a vertex and a leg. 
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diagram 5F — 5G 

3. (a) 

e2 

5. (a) - — p-p2T(p,p2,Mi^, M^rjlnm^ 

5. (b) -^p-piT{p,pi,MK,M^)]nm'^ 

e2 

5. (c) —p ■ piTip,pi, MK,mi)\nm^ 

5. (d) -— pi •p2T"(pi, -p2,M^,il4) Inm^ 
e2 

5. (e) —pi-piT{pi,-pi,M^,mi)liim'^^ 
e2 

^- (^) -—^P2-Pir{p2,-pi,MT„mi)]nm'i 



Table 1: Infrared divergent part of the corrected / and g form factors due 
to virtual photon corrections. The contribution from diagram 3. (a) comes 
from wave function renormalization of external charged particles, K~^, tt"*", 
7r~, and 
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/+ /+ 

(b) 



Figure 7: Virtual photons in Feynman diagrams for the decay in the 
charged channel. One takes the pure vertex mediated by a K'^ pole and then 
attaches a photon between two legs. 
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